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Abstract 

The disorder function formalism [Gunaratne et.al., Phys. Rev. E, 57, 5146 (1998)] is used to 
show that pattern relaxation in an experiment on a vibrated layer of brass beads occurs in three 
distinct stages. During stage I, all lengthscales associated with moments of the disorder grow 
at a single universal rate, given by L(t) ~ t 5 . In stage II, pattern evolution is non-universal 
and includes a range of growth indices. Relaxation in the final stage is characterized by a single, 
non-universal index. We use analysis of patterns from the Swift-Hohenberg equation to argue 
that mechanisms that underlie the observed pattern evolution are linear spatio-temporal dynamics 
(stage I), non- linear saturation (stage II), and stochasticity (stage III). 
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I. INTRODUCTION 



Analysis of structures far from equilibrium and their evolution has necessitated the devel- 
opment and implementation of several new techniques jlj]. The underlyingphysical systems 
are dissipative and typically exhibit sensitive dependence on initial states |2(; i.e., differences 
between nearly identical patterns increase in time. Further, although structures formed from 
distinct initial states differ in detail, they appear to share features that depend only on con- 
trol parameters of the physical system. These observations led Prigogine and collaborators 
to propose the use of irreversible statistical analysis using ensembles of solutions of a sys- 
tem [3j, y, la, |y] . However, it is a non-trivial task to identify measures appropriate to analyze 
a given class of structures. 

We consider properties of large aspect ratio patterns that are formed via spontaneous 
symmetry breaking from a noisy, but otherwise structureless initial state. Patterns are gen- 
erated when an external control parameter crosses a threshold, initiating a local bifurcation 
to a periodic array; these arrays are typically, striped, square, triangular or hexagonal jlj]. 
However, due to translational and rotational invariances of the underlying system, the lo- 
cation and orientation of these local patches are not synchronized over the entire domain. 
As they grow, these patches compete with each other to produce complex patterns seen in 
experimental and model systems 1 1 31 ] , see Figures ^ and El These 

structures, whose details are sensitively dependent on the initial state, can be considered to 
belong to the ensemble of patterns. 

Assume that such a spatio-temporal pattern can be represented by a scalar field u(x,t). 
Several functionals of w(x, t) have been proposed as measures to be used in statistical de- 
scriptions. Most of these are modeled on related quantities from equilibrium statistical 
mechanics. The width of the structure factor S(t) jla . 14, ^j (i.e., the width of the peak 
in the azimuthal average of (w(k, t)u(— k, t))) is one of the more frequently studied exam- 
ples. It has been shown that the evolution of a pattern from a noisy initial state consists of 
two distinct stages [id ]: S(t) ~ is obeyed until the peak amplitude of the field u(x, t) 
saturates, beyond which time the pattern coarsens and the decay becomes slower. For the 
Swift-Hohenberg equation, Elder et al. [l^, Cross and Meiron 15], and Hou et al. [3] found 
that during this second stage S(t) decreased as t~ a20 when the dynamics was noise free and 
at a faster rate t~ - 25 otherwise. Additional quantities have been introduced to describe 
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the second relaxation region. Hou et al. |3| measured cumulative domain wall length as 
a function of time. For noise-free (resp. noisy) dynamics, it was shown to decay as t~ a25 
(resp. t~ a3 °). Cross et al. .15] showed that a stripe orientation correlation field decayed as 

r 0.24 

Notice that in the second relaxation stage, distinct measures (which probe multiple fea- 
tures of a pattern) decay at different rates. A systematic way to study such variations is to 
define a field to quantify the "local" disorder, and to use multiple moments of this field to 
construct a family of characteristics to describe the dynamics of various lengthscales asso- 
ciated with the pattern [17| . Such a class appropriate for labyrinthine patterns, referred to 



as the disorder function, has been introduced Q • It is defined as 

'(2-[3)fd 2 x\(A + k 2 )u(x,t)\^ 



8(M 



k 2 ?, (1) 



(jd 2 x) (\u(x,tW 

where k is the wavevector associated with the local planform, (|iz(x, t)\) denotes the mean 
of |w(x, t)\, J d 2 x is the area of the pattern, and S(j3,t) (0 < (3 < 2) has been normalized 
to be scale invariant. In this definition, the field used to quantify the local disorder is 
the magnitude |(A + k 2 )u{x, t)\ of the Helmholtz operator. The integrand is non-zero at 
locations where the local wave-number deviates from its global mean k, and those where the 
planforms are curved. (The latter can be seen by considering a region of stripes with a radius 
of curvature r; i.e., w(x, t) ~ exp(ikr). There |(A + k 2 )u(x, t)\ ~ k/r.) Distinct "moments" 
(3 can be used to emphasize particular features of a pattern; for example, large values 
of (3 provide a higher weight to defects |3, I20I1 . Methods to calculate S(/3,t) for (noisy) 
experimental patterns have been introduced [2JJ. In this paper, we will describe spatio- 
temporal dynamics using 5'(f3,t) = 5 ((3, t) 1 ^. As can be easily seen, for each (3, 5'((3,t) is 
the ratio of a wave number associated with the moment (3 of the Helmholtz operator and k. 

The disorder function analysis of the evolution of patterns under model spatio-temporal 
dynamics from an initially noisy state reveals the presence of two stages Q], consistent with 
the analysis of S(t) [13]. It is found that during the first stage spatio-temporal dynamics is 
essentially linear, and 5'(f3,t) ~ t~~ ai , where oj = | is a universal exponent [19]. Thus, the 
wave number associated with each moment f3 reduces at a rate t~ UI , and the corresponding 
lengthscale grows as t ai . The second stage, domain-coarsen ing, exhibits more complex 



behavior. The relaxation is described by 5'((3,t) ~ t CTlI ^' t > 20(. Unlike 07, au((3,t) is 



model dependent; i.e., the corresponding relaxation is not universal, and wave numbers 



associated with multiple moments decay with rates that depend on (3. 

These conclusions were recently validated in a set of experiments on pattern formation 
in a vibrated layer of brass beads j^J. However, the experiments showed, in addition, 
that for very long times the patterns reach a third stage; here, the indices crju(/3,t) are 
independent of /3, but their common value is not universal. The work reported in this paper 
shows that the presence of this final stage in systems far from the onset of patterns is one 
observable consequence of stochasticity in the underlying spatio-temporal dynamics. The 
results com plem ent recent demonstrations on effects of noise in systems close to the onset 
of patterns 
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II. MODEL SYSTEM AND ANALYSIS 

The model system used in our studies is an extension of the Swift-Hohenberg equation Q 
which describes the dynamics of a scalar field w(x, t) via 

du 



^ = £>( e - (A + klf)u -u 3 - v(Vuf + r/(x, t). (2) 

D sets the rate of diffusion, e is the distance from pattern onset, v is the strength of a 
non- variational term 3|, and ko is the wave number of the underlying planform. rj is a 
stochastic term such that (?7(x, t)i](x.', £')) = F6(x — x.')6(t — t'), F being the strength of the 
noise; rj 7^ and rj = represent stochastic and noise-free spatio-temporal dynamics. For 
suitable control parameters, random initial states (whose amplitude is chosen to be much 
smaller than the saturated peak amplitude of the field) evolve to patterns under Eqn. (J2J)- 
Several snapshots from such an evolution under noise-free dynamics are shown in Figure Q 
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The behavior of 5(1.0, t) corresponding to this spatio-temporal dynamics is shown in 
Figure|21 Statistical variation between values of 6(1.0, t) in nominally identical patterns-i.e., 
those with same control parameters, but different initial states- is neglig ible Q. 6(1.0, t) 
exhibits an initial rapid decay followed by a much slower relaxation. During stage I, (\u(x,t)\) 
is significantly smaller than its value at large times; hence, spatio-temporal dynamics can 
be approximated by the linear partial differential equation 

^ = (e - (A + klf)u. (3) 



5 



The Fourier transform of w(x, t) satisfies 

u(k, t) = u(k, 0) exp{et - (jfcg - k 2 ) 2 t). (4) 

For a given noisy initial state, this solution can be used to evaluate integrals in Eqn. (JTJ). It 
is found that (\u(x,t)\) ~ e et t _1//4 and (for values of t that are not very small) / d 2 x\(A + 



k*)u(x,t)\P ~ e^H-^l 4 25]. Consequently, 8({3,t) ~ r# 2 and cr(/3, t) = a l = 1/2. 

Results shown in Figure El can be used to compute a(/3,t). At a given time t, we use 
the five nearest points (t and the closest two points on both sides) to calculate the slope of 
the curves log 5 (ft, t) vs. logt using a least squares fit. As shown in Figure 01 this analysis 
confirms that 0/(/3, t) = 1/2 during stage I. 

The onset of stage II coincides with the saturation of the peak value of the field u(x, t) and 
beyond this time the relaxation of 5(1.0, t) is slower. Further, as seen from Figure E3 (r(/3,t) 
varies both with t and (3. The variation with (3 implies that the wave numbers associated 
with 8'(/3,t) decay at rates which depend on j3\ i.e., lengthscales associated with multiple 
features of these patterns grow at different rates (^(J. The spatio-temporal dynamics shown 
in Figure retains this feature even for very large times, when the pattern consists of a few 
large domains. 

Typically, a(/3, t) is a monotonically decreasing function of j3 during stage II; thus, the 
mean distance between defects (/3 — > 2 jl^) grows at a slower rate than (for example) the 
mean curvature or domain size. Exceptions occur during the removal of defect (s) or small 
domains from a pattern. These events can be identified by a sudden rapid decrease in 8(/3, t), 
an example of which is seen near T4 in Figure El Figure 0] shows the pattern prior to and 
following the absorption of a small domain into a larger one. The behaviors of 8(1.0, t) and 
a(/3,t) during this metamorphosis is shown in Figure 02a) and (b). It can be seen that 
cr(/3, t) which is monotonic before and after the change takes on a more complex form during 
the transition. This method can be used to identify domain disappearance even when it is 
difficult to recognize the change visually (as often happens in patterns with a large number 
of domains). 
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III. RESULTS FROM THE EXPERIMENT 



The pattern forming experiments were conducted with 0.165 mm bronze spheres con 
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tained in a vertically oscillated circular container with a diameter of 140 mm 
layer is four particle diameters deep, and the cell is evacuated to 4 Pa so that hydrodynamic 
interaction between the grains and surrounding gas is negligible. The control parameters 
are the frequency / of the sinusoidal oscillations and the peak acceleration of the container 
relative to gravity, T = (27rf) 2 A 2 /g, where A is the amplitude of the oscillation and g is the 
gravitational acceleration. As / and F are changed, a variety of temporally subharmonic 
patterns including locally square, striped, or hexagonal patterns are observed jllj]. The 
textures analyzed consisted of patterns with square planforms j^ . 

The granular surface is visualized by illuminating the cell using a ring of LEDs surround- 
ing the cell and is strobed at the drive frequency of the container. The light is incident at 
low angles and the scattering intensity is a nonlinear function of the height of the layer; 
scattering from peaks (valleys) creates bright (dark) regions. This intensity field is used to 
represent the structures j^J. 

Figure shows four snapshots during the evolution of an experimental pattern, and 
Figure [7| shows the corresponding behavior of 5(1.0, t). Stages I and II, analogous to those 
seen in patterns of the model system, can be observed. It was confirmed that the onset 
of stage II corresponded to the saturation of the peak amplitude j^j]. Using the leftmost 
interval shown in Figure the evolution during stage I was shown to be consistent with 
<j(j3, t) = 1/2, although there were insufficient data points for the result to be conclusive J^. 
Beyond this stage, the relaxation was much slower, and a(/3,t) depends on /3 jscj, see 
Figure |H1 However, it was observed that at much longer times (e.g., during the rightmost 
interval in Figure EJ) the patterns reach a stage where a(/3,t) is independent of /3, and its 
value is not the universal 1/2 of the first stage, see Figure |H1 



IV. SPATIO-TEMPORAL DYNAMICS UNDER STOCHASTIC EQUATIONS 

Behavior analogous to stage III in experimental patterns is not observed in the spatio- 
temporal dynamics such as that shown in Figure even at very large times when the 
pattern consists of a few domains. Clearly, one or more ingredients that are necessary for 
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the description is missing from the model we have used. Candidates include non-variational 
terms, finite size effects and stochastic terms. In particular, only finite subdomains of the 
experimental system are analyzed (in contrast the model system has periodic boundary 
conditions), and stoehastic.ty has been shown to play a to.e in the dynam.es neat the onset 
of patterns 23] . 

Integration of the model system shows that neither the addition of non-variational terms 
(v ^ 0) to the spatio-temporal dynamics nor conducting the analysis on smaller sub-domains 
give a transition away from stage II, even at very large times. On the other hand, we find 
that the addition of stochastic terms does lead to such a transition, and that the new stage 
has characteristics described in the last section. 

Figure El shows four snapshots during the development of such a structure under the 
Swift-Hohenberg equation with stochastic terms (i.e., F 7^ 0). The behavior of 5(1.0, t) 
is qualitatively similar to Figure 121 The principal difference that can be gleaned from the 
analysis of a single moment 5((3,t) is an increase in the decay rate, particularly in stage 
II. a(/3,t) exhibits features that are similar to the corresponding noise-free dynamics; it is 
independent of (3 in stage I, but is a function of (3 following saturation of the field w(x, t). 
However, at very large times, a(/3,t) is once again seen to be independent of (3, taking a 
value that depends on control parameters, see Figure [TU1 It is also found that the onset of 
stage III advances with increasing intensity of the stochastic terms. For sufficiently large 
values of F the dynamics can move directly from stage I to III. 

The implication of the j3 independence of <Jni is that all lengthscales associated with the 
pattern grow at the same rate during stage III. This assertion can be tested by studying 
the relaxation of the width of the structure factor and the defect density. For the analysis, 
we used averages for four nominally identical runs with parameters given in Figure EJ The 
mean decay rate of the width of the structure factor at time T 3 of Figure is found to be 
0.119 ± 0.014, in agreement with <tui(T 3 ). Secondly, the mean number of defects at T 3 (XX 
HOW WAS THIS CALCULATED?) is found to decrease like T x , where x = -0.201 ±0.019. 
The associated lengthscale (i.e., the mean distance between defects) grows with an index 
|%, once again consistent with the value of amiT^). 
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V. DISCUSSION AND CONCLUSIONS 



We set out to compare and contrast pattern dynamics in distinct systems. It is im- 
portant to note that such a task can only be achieved using "configuration independent" 
characteristics; i.e., measures that depend on control parameters of the underlying system, 
and not on the (random) initial states. For complex patterns with local planforms such as 



stripes, squares or hexa gon s (e.g., Figures 



the structure factor 
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lq . defect density 



ElandlHl), such measures include 
lq . length of domain boundaries 



;he width of 
I, and the 



disorder function |l8l |29|. 

Analysis of patterns generated in the Swift-Hohenberg equation and on a vibrated layer 
of brass beads show that during early times all components of the disorder function relax as 
£-0.5 ^y e ghow-gc! that the index can be obtained by noting that during this stage the intensity 
of the field representing the pattern is sufficiently small and hence the spatio-temporal 
dynamics can be considered to be linear. The saturation of the field is followed, in each 
case, by a second stage where the relaxation is non-universal and the growth indices a(/3,t) 
are (3 dependent; the latter implies that the growth of distinct features of patterns occur at 
different rates. Experimental patterns exhibit a third stage at times very long compared to 
the vibrational period of the container where a(/3,t) is once again independent of j3. It 
is not possible to differentiate it from stage II by analyzing the behavior of a single index; 
rather, it requires the use of the entire family of indices a(/3,t). Unlike in stage I, the value 
of a(P, t) is not universal. Stage III is not present in spatio-temporal dynamics under the 
non-stochastic Swift-Hohenberg equation, even when non- variational terms and finite-size 
effects are included. Our work suggests that the ingredient required to initiate stage III is 
stochasticity of the underlying spatio-temporal dynamics. 

It should be noted that 5(/3,t) was constructed to quantify disorder in patterns that 
consist of local planforms which satisfy the Helmholtz equation [29]. Structures formed in 
other systems, such as in the X-Y model pOI ]. in spinodal decomposition [^jI, or in epitaxial 
growth [32] do not have such local planforms. Measures to conduct irreversible statistical 
analyses of such structures remain to be identified. 

The authors would like to acknowledge discussions with M. Shattuck. This research is 
supported by grants from the National Science Foundation (DJK, GHG), the R. A.. Welch 
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FIG. 1: Four snapshots of the spatio-temporal dynamics of a random initial state (|it(x, t)\ < 0.001) 
under Eqn. E] with D = 0.01, e = 0.1, k = 3, and F = 0. Times T x = 124 and T 2 = 924 are in 
stage I, and T 3 = 42, 528 and T 4 = 939, 878 are in stage II. 
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FIG. 2: The behavior of 6(1.0, T) for the spatio-temporal dynamics shown in Figure^ The presence 
of two distinct stages of relaxation can clearly be observed. 
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FIG. 3: The behavior of the function <r(/3, T). During the domain forming stage cr(/3, T) = 1/2, 
while it is /? dependent in the domain coarsening stage. As in the example shown, ct(/3, T) also 
varies with time during the second stage. 




FIG. 4: Sharp reductions in 5{(5,T) (similar to that near T4 in Figure ^) correspond to the disap- 
pearance of domains or defects. The images (a) before, and (b) after such an event. It is found 
that <j{(5, T) fails to be monotonically decreasing during and immediately following such events. 
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FIG. 5: The behavior of 5((3, t) and cr(/3, t) during the removal of a small domain. Notice that, while 
a((3, t) decreases monotonically before the event and following the relaxations of the neighborhood, 
it exhibits more complex behavior during the event. 
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45 100 

FIG. 6: the time evolution of an oscillated layer of grains following a rapid change in peak accelera- 
tion of the container from 2.2g to 3.0g. Below the onset of patterns (which occurs at approximately 
2.5<7) 5 the layer is featureless but noisy while above onset defect-free square patterns are stable. 
The numeral below each snapshot indicates the number of layer-plate collisions after the last layer 
takeoff before the acceleration was changed. The oscillation frequency in this experiment is / = 27 
Hz. 
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FIG. 7: The time evolution of 5(1.0) for the sequence shown in Figure H3 The abscissa is in units 
of the number of container oscillations from the last layer takeoff take before T was changed. The 
dashed lines indicate the time that snapshots in Figure|H]were taken. Two distinct stages of pattern 
evolution are clearly visible; the initial rapid decay is consistent with i -0,5 , while the slower decay 
scales on average like i -0 ' 12 . However, it is difficult to differentiate stages II and III from the 
dynamics of a single moment of the disorder function. 
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FIG. 8: a(P,T) calculated for patterns shown in Figure H3 Results from stages I, II, and III are 
shown in circles, diamonds and squares respectively. <r(/3, t) for the three stages were calculated us- 
ing intervals shown in Figure |H] (3-5 oscillations, 40-50 oscillations, and 85-105 oscillations following 
the initiation of the experiment). Now stages II and III can clearly be delineated. 
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FIG. 9: Four snapshots of the stochastic spatio-temporal dynamics of a random initial state 
(\u(x,t)\ < 0.001) under Eqn. |U with e = 0.01, k = 3, and F = 0.0001. The T x = 924 Is in 
stage I, T2 = 45,528 in stage II, and the last two T3 = 546,023 and T4 = 1,630,417 are in stage 
III. Note that the pattern continues to evolve, albeit slowly, in stage III. 
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FIG. 10: The behavior of a((3, T) for stochastic spatio-temporal dynamics. Observe the presence 
of three distinct stages, as in experimental patterns. 
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